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Modifying and extending recent ideas i, a theoretical framework to describe dissipation processes 
in the surfaces of vibrating micro and nanoelectromechanical devices (MEMS-NEMS), thought to 
be the main source of friction at low temperatures, is presented. Quality factors as well as frequency 
shifts of flexural and torsional modes in doubly-clamped beams and cantilevers are given, showing 
the scaling with dimensions, temperature and other relevant parameters of these systems. Full 
agreement with experimental observations is not obtained, leading to a discussion of limitations and 
possible modifications of the scheme to reach quantitative fitting to experiments. For NEMS covered 
with metallic electrodes the friction due to electrostatic interaction between the flowing electrons 
and static charges in the device and substrate is also studied. 

PACS numbers; 03.65.Yz, 62.40.+i, 85.85.+j 



I. INTRODUCTION 

The successful race for miniaturization of semiconduc- 
tor technologies — manifests itself spectacularly in the 
form of nanoelectromechanical systems (NEMS) — lii^i^, 
machines in the micron and submicron scale whose me- 
chanical motion, integrated into electrical circuits, has a 
wealth of technological applications, including control of 
currents at the single-electron level — , single-spin detec- 
tion sub-attonewton force detection — , mass sensing of 
individual molecules high-precision thermometry — 
or in-vitro single-molecule biomolecular recognition — . 

These mechanical elements (like cantilevers or beams, 
see fig. [J) are also a focus of attention and intensive 
research as experiments are approaching the quantum 
where manifestations of quantized me- 
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regime 

chanical motion of a macroscopic degree of freedom like 
their center of mass should become apparent. Several 
schemes to prepare the mechanical oscillator in a non- 
classical state and observe clear signatures of its quantum 
behavior have been recently suggested 16^17,18^19,20^21 ^ 

A key figure of merit of the mechanical oscillation is 
its quality factor Q = cj/Aw, where Auj is the measured 
linewidth of the corresponding vibrational eigenmode of 
frequency uj. To reach the quantum regime, as well as for 
most practical applications, where the measured shifts of 
the resonant frequency u constitute the detectors princi- 
ple, a very high Q is compulsory. Imperfections and the 
environment surrounding the oscillator result in both a 
finite linewidth Au> and a frequency shift Slo with respect 
to the ideal case. Therefore several works have been de- 
voted to the analysis of the different sources of dissipation 
present in MEMS and NEMS ^^^^^^^28^^^^^ 
trying to determine the dominant damping mechanisms 
and ways to minimize them. Among the different mech- 
anisms affecting semiconductor-based NEMS the most 
important and difficult to avoid are i) clamping losses 



FIG. 1: (Color online) Sketch of the systems considered in 
the text. The inset shows a doubly clamped beam, while the 
main figure shows a cantilever characterized by its dimensions, 
width (w), thickness (t), and length, (L), where w ^ t <^ L. 
The height above the substrate is h. A schematic view of the 
surface is given, highlighting imperfections like roughness and 
adsorbates, which dominate dissipation at low temperatures. 



^iM, through the transfer of energy from the resonator 
mode to acoustic modes at the contacts and beyond to 
the substrate, ii) thermoelastic damping ^^i'^^'^^ and iii) 
friction processes taking place at the surfaces 38,39,4o _ p^^ 
low temperatures and for decreasing sizes the prevailing 
mechanism is the last one — , as indicated by the linear 
decrease of the quality factor of flexural modes with de- 
creasing size (see fig. [5]), or the sharp increase of Q when 
the resonator is annealed — i^i. Excitation of adsorbed 
molecules, movement of lattice defects or configurational 
rearrangements absorb irreversibly energy from the ex- 
cited eigenmode and redistribute it among the rest of 
degrees of freedom of the system. 

A theoretical quantitatively accurate description of 
surface dissipation proves therefore challenging, as many 
different dynamical processes and actors come into play, 
some of whom are not yet well characterized, so simplifi- 
cations need to be made to provide a unified framework 
for all of them. In such a scheme was given, based 
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FIG. 2; (Color online) From ref. ^. Evolution of reported 
quality factors in monocrystalline mechanical resonators with 
size, sho-wing a decrease with linear dimension, i.e., with in- 
creasing surface-to- volume ratio, indicating a dominant role 
of surface-related losses. 



on the foUo'wing considerations: i) experimental observa- 
tions indicate that surfaces of other'wise monocrystalline 
resonators acquire a certain degree of roughness, im- 
purities and disorder, resembling an amorphous struc- 
ture — . ii)In amorphous solids the damping of acous- 
tic 'waves at low temperatures is successfully explained 
by the Standard Tunneling Model isdid^d^, which cou- 
ples the acoustic phonons to a set of Two-Level Sys- 
tems (TLSs) representing the lo'w-energy spectrum of all 
the degrees of freedom (DoF) able to exchange energy 
with the strain field associated to the vibration. These 
DoF correspond to impurities or clusters of atoms within 
the structure which have in their configurational space 
two energy minima separated by an energy barrier (sim- 
ilarly to the dextro/levo configurations of the ammonia 
molecule), modeled as a DoF tunneling between two po- 
tential wells. At low temperatures only the two low- 
est eigenstates have to be considered, characterized by 
the bias Ag between the wells and the tunneling rate 
Ag through the barrier. The Standard Tunneling Model 
specifies the properties of the set of TLSs in terms of a 
probability distribution ^(Aq, Ag) which can be inferred 
from general considerations — lii and is supported by ex- 
periments — . It has to be noted, however, that below a 
certain temperature the model breaks down due to the 
increasing role played by interactions among the TLSs, 
not included iliiS. 

In — a description of the attenuation of vibrations in 
nanoresonators due to their amorphous-like surfaces in 
terms of an adequate adaptation of the Standard Tun- 
neling Model was intended to be given. An estimate for 
Q^^{T) was provided, reproducing correctly the weak 
~ T^/^ temperature dependence as well as the order of 
magnitude observed in recent experiments ^i^. How- 
ever, the way the Standard Tunneling Model was adapted 
was not fully correct, and a revision of the results for 



Q i(r) obtained therein was mandatory. 

In this work we will i) Discuss in detail some impor- 
tant issues hardly mentioned in ^ , ii) Modify some points 
to obtain a theory fully consistent with the Standard 
Tunneling Model and which includes all possible dissi- 
pative processes due to the presence of TLSs, iii) Extend 
the results and give expressions for the quality factor of 
cantilevers as well as damping of torsional modes, and 
frequency shifts associated, iv) Compare with available 
experimental data, discussing the validity of the model, 
range of applicability, aspects to be modified/included 
to reach an accurate quantitative fit to experiments, v) 
Study the dissipative effects associated to the presence 
of metallic electrodes frequently deposited on top of the 
resonators. 

Section II starts with a brief summary of the model 
described in ^, discussing the approximations involved, 
while the different dissipative mechanisms due to the 
presence of TLSs are presented in III. Section IV ana- 
lyzes and compares the two main mechanisms, namely 
relaxational processes associated to asymmetric (biased) 
TLSs, and non-resonant damping of symmetric TLSs. 
The extension of the results to cantilevers and torsional 
modes is given in Section V. A brief account of frequency 
shifts' expressions is found in Section VI. In Section VII 
comparison with experiments and discussion of the ap- 
plicability, validity and further extensions of the model 
is made. Section VIII discusses dissipation effects due to 
the existence of metallic leads coupled to the resonator, 
which can be coupled to the electrostatic potential in- 
duced by trapped charges in the device. Under some cir- 
cumstances, these effects cannot be neglected. Finally, 
the main conclusions of our work are given in Section 
IX. Details of some calculations are provided in the Ap- 
pendices. 

We will not analyze here dissipation due to clamping 
and thermoelastic losses, which may dominate dissipa- 
tion in the case of very short beams and strong driving, 
respectively. We will also not consider direct momentum 
exchange processes between the carriers in the metallic 
circuit and the vibrating system — . 



II. SURFACE FRICTION MODELED BY TLSs 



Hamiltonian 



We will consider a rod of length L, width w and thick- 
ness i, see fig. [I], and use units such that Ti ~ I = ks- 
As described in — , the vibrating resonator with imper- 
fect surfaces is represented by its vibrational eigenmodes 
coupled to a collection of non-interacting TLSs, assuming 
that the main effect of the strain caused by the phonons 
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is to modify the bias Aq of the TLSs 



50. 



(1) 



k,3 



The index j represents the three kinds of modes present 
in a thin beam geometry — : flexural (bending), tor- 
sional and compression modes. These modes will be 
present for wavelengths A > while for shorter ones 
the system is effectively 3D, with the corresponding 3D 
modes. The main effect of these high energy 3D modes is 
to renormalize the tunneling amplitude, so we will take 
that renormalized value as our starting point Aq and 
forget in the following about the 3D modes. The sum 
over TLSs is characterized by the probability distribu- 
tion P(A§, A 
explained in 

tivity of Ag to the properties of the energy barrier of the 
two-well potential giving rise at low T to the TLS de- 
scription of the system (resulting in a 1/ Ag dependence 
of -P(Ag, Ag)), and ii) the characteristic energy scale of 



Q, - Po/A§ . i^ii This result follows, as 
— 1^', just due to i) the exponential sensi- 



the distribution of asymmetries Aq, much bigger than 
IK, which is the temperature at which experiments are 
performed, and which thus fixes the scale of the bias of 
a TLS if it is to contribute significantly to dissipation, 
A§ < IK (therefore the relevant TLSs have values of Ag 
lying in a very narrow energy range around Aq = as 
compared to the variance of their probability distribu- 
tion, allowing us to consider that P(Aq,Aq) to a first 
approximation does not depend on Ag). Unphysical di- 
vergencies do not appear, as Ag > Ami„, with A^m fixed 
by the typical timescale of the experiment, given by the 
time needed to obtain a spectrum around the resonance 
frequency of the excited vibrational eigenmode of the res- 
onator, and e — (Ag)^ -|- (Ag)^ < Smaxi estimated to 
be of the order of 5 K — . For typical amorphous insu- 
lators Po 1044j-im-3. 

To see to which low energy modes the TLSs are 
more coupled, inducing a more effective dissipation 
at low temperatures, the spectral function J(a;,j) = 

X]fe l-^fe,i I ^{'^ ~ ^k,j) characterizing the evolution of the 
strength of the coupling for each type of mode can be 
computed. Due to their nonlinear dispersion relation, 
Lu = ^/WJ{ptw)P, with / = t^w/12, E the Young mod- 
ulus and p the mass density, flexural modes show a sub- 
ohmic behavior Jflex(w) = a\,JTJ^^fuj^ with 



= 0.3 



7^ {\ + v){\-2v) fp 



1/4 



(2) 



where 7 ~ 5 eV is a coupling constant appearing in 
Xkj, V is Poisson's ratio and Wco — \J EI j {ptw){2'K /tY 
is the high energy cut-off of the bending modes. A de- 
tailed derivation of Jflcx(^) is given in Appendix [XI Even 
though the length L of our system is finite, and thus the 
vibrational spectrum discrete, a continuum approxima- 



tion like this one will hold if kT ^ fiufund^ ^fund being 
the frequency of the fundamental mode. 

The bending modes prevail over the other, ohmic- 
like, modes as a dissipative channel at low ener- 
gies, thanks to their weaker Jflox(^) ~ uj^/"^ depen- 
dence. One may ask at what frequency do the tor- 
sional and compression modes begin to play a signif- 
icant role, and a rough way to estimate it is to see 
at what frequency do the corresponding spectral func- 
tions have the same value, Jflex('^*) = '^comp,tors('i'*)- 
Using the expressions ini, namely Jcomp,tors(^) — 
a^.tuj, with ac = {jA^/Aof{2Tr'^ptw)-^E/p)-^^^ and 
at = C(7Ag/AQ)2(87rVtwp/)"^(p//C)3/2, the results 
are cu* - 30(1 + 1^)^(1 - 2iy)^{E/pY/y[t{3 ~ 5i^)^] for 
the case of compression modes and uj* ^ 300(1 — 
2zy)2(£;/p)i/2/[i(3 _ + „)] for the torsional. Com- 

paring these frequencies to the one of the onset of 3D 
behavior, lUco, they are similar, justifying a simplified 
model where only fiexural modes are considered. 



B. TLSs dynamics 

The interaction between the bending modes and the 
TLSs affects both of them. When a single mode is ex- 
ternally excited, as is done in experiments, the coupling 
to the TLSs will cause an irreversible energy flow, from 
this mode to the rest of the modes through the TLSs, as 
depicted in fig.([3K). The dynamics of the TLSs in pres- 
ence of the vibrational bath determines the efficiency of 
the energy flow and thus the quality factor of the ex- 
cited mode. Taking a given TLS plus the phonons, its 
dynamics is characterized by the Fourier transform of the 
correlator {az{t)(Jz{0)), the spectral function A{llj), which 
at T = reads 



A(a;)^^|(0|a, 



■d{Lu- 



(3) 



where \n) is an excited state of the total system TLS plus 
vibrations. In — an analysis of A(uj) was made for the 
case of a symmetric TLS (Ag = 0), concluding that i)If 



Ag <C a^Wco the tunneling amplitude is basically sup- 



pressed and the TLS does not participate in dissipative 
processes. For reasonable system dimensions the cou- 
pling constant is very small, at 1, so this effect can be 
ignored, ii)Around the resonance at lu ~ Ag a broaden- 
ing appears, r(Ag) which, for r(Ag) <C Ag is given by 
the Fermi Golden Rule result r(Ag) — X'oahsJiOco \/ Ag 
(for T = K, at T > eq.© appfies), iu)The coupling to 
phonons of all energies provides the "dressed" TLS with 
tails far from resonance, A{lo) oc af,^/cj^^^/(Ag)^ for 
w < Ag and A(w) cx af,^/ZI^(Ag)2w-7/2, for w > Ag, 
see left side of fig.(13b), iv) The main effect of the asym- 
metry Ag is to suppress the TLSs dynamics, so the TLSs 
playing an active role in dissipation satisfy Ag > |Ag|. 
Finally, for those low-energy overdamped TLSs such 



that r(Ag) > Ag 
time, with T(Ag) : 



an incoherent decay exp(— t/r) with 
= r(Ag)~^ was assumed. 
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FIG. 3: (Color online)a) Schematic representation of the ir- 
reversible flow of energy from the externally excited mode to 
the TLS ensemble, and from the ensemble to the rest of vibra- 
tional modes. This process can be viewed as a flow of energy 
from the excited mode to an ensemble of "dressed" TLSs, with 
their dynamics modifled by the presence of the vibrational 
modes, b) Left: Spectral function A(yj) of a single dressed 
TLS, weakly damped (F < Ag). A peak around a; = arises 
if the system is biased, corresponding to the relaxational dis- 
sipation mechanism. Right: Total spectral function Atoti}^) 
of the ensemble of dressed TLSs. 



Two points that ref. — missed, and now will be con- 
sidered, are the following: First, the effect of the asym- 
metry is not simply the one stated. An underdamped 
biased TLS develops an additional peak around uj = 0, 
as shown in fig.llSb), which corresponds to the relaxation 
mechanism that dominates dissipation of acoustic waves 
in amorphous solids — , see eq.®. In — the relaxation 
mechanism associated to eq.® was misinterpreted to 
correspond to friction due to overdamped TLSs, while 
in this work it will be, consistently with the Standard 
Tunneling Model, linked to the presence of biased under- 
damped TLSs. It will be described in more detail in the 
next section, and taken into account in the computation 
of the total dissipation. 



Second point missed by ref. ^: one can estimate, using 
the probability distribution P(Ag,AQ), the total num- 
ber of overdamped TLSs in the volume fraction of the 
resonator presenting amorphous features, Vamorph- With 
Vamorph ~ Vtot/10, and usiug T{E,T) - T{E,T - 
0) coth[£'/2r] « 2TT{E, T = 0)/E, the number of over- 
damped TLSs, Ai5 < r(Ag,r) ^ Ag < [30a67^I^r]2/3, 
is « PotwL [SOab^/uJcoT] ^/"^ , which for typical res- 
onator sizes L ~ 1/xm, t,w ^ 0.1/im is less than one 
for T < IK. Therefore, unless the resonator is bigger 
and/or Pq too, this contribution to dissipation can be 
safely neglected, as it will be done in the following. 



III. DISSIPATIVE MECHANISMS 

From the features of the spectral function Atot{(^) — 
X]a=" ^(^O'^O''^) of the ensemble of dressed TLSs, 
one can classify into three kinds the dissipative mecha- 
nisms affecting an externally excited mode lvq: 

A. Resonant dissipation 

Those TLSs with their unperturbed excitation ener- 
gies close to Wo will resonate with the mode, exchanging 
energy quanta, with a rate proportional to the mode's 
phonon population riuj^ to first order. For usual excita- 
tion amplitudes ^ lA the vibrational mode is so pop- 
ulated (as compared with the thermal population) that 
the resonant TLSs become saturated, and their contri- 
bution to the transverse (flexural) wave attenuation be- 
comes negligible, proportional to 



'1/2 



52 



B. Dissipation of symmetric non-resonant TLSs 

We show first that a correct description of the dissi- 
pation due to weakly damped TLSs is given by the ap- 
proximate Hamiltonian used in ^, starting from eq.([T]), 
if one considers apart in some way the presence of the 
relaxational peak at w = due to a finite bias Aq ^ 0. 

We focus the attention on a given TLS plus the vibra- 
tions (spin-boson model), H = Agcr^: + AqCT^ -|- Hint + 



Hint — (^z ^ A 



fc2 



/LUk 



:(4 



(4) 



with A defined in eq. (|A4p . In — a change of basis to the 
unperturbed eigenstates of the TLS was performed, ob- 
taining H= sa, + [{A-/e)a, + (Ag/e)a,] A^(al + 
ak) + Hyihr- Then the term {AQ/e)az was ignored (re- 
member £ = ■\/(Ag)^ + (Ag)2 ) , arguing that the key 
role in dissipation is played by fairly symmetrical TLSs, 
\Aq\ ^ e. This allowed to simplify the spin-boson Hamil- 
tonian to the one of a symmetric TLS, much easier to 
analyze, with tunneling amplitude e instead of Ag and 
coupling A(Ag/e) instead of A. The consistency was kept 
by restricting the sums over the TLS ensemble to those 
such that |Ag| < Ag. 

One can check this consistency going back to the 
original basis, where the approximation of the Hamil- 



[(Ag/£) 



tonian reads H « Agcr^; 

{A^A^/e^)ax]J2k^ (al + a/c ) + H^ibr ■ Restricting the 
application of this Hamiltonian to those TLSs such that 
I Agl < Ag seems to be a fairly good approximation, but a 
price has been paid, namely the spectral weight at = 
due to the bias has been lost in this effectively symmet- 
ric spin-boson approximation. This weight cannot be ig- 
nored, and it has to be added as a different mechanism 
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(the relaxational mechanism), what will be done in the 
next subsection. Once this issue has been taken care 
off, all the dissipative processes due to the presence of 
TLSs are correctly included in this framework, and we 
can proceed describing non-relaxational friction due to 
non-resonant, weakly damped, weakly biased TLSs. 

As mentioned before, the coupled system TLSs + vi- 
brations can be viewed, taking the coupling as a pertur- 
bation, from the point of view of the excited mode ojo as 
a set of TLSs with a modified absorption spectrum. The 
TLSs, dressed perturbatively by the modes, are entities 
capable of absorbing and emitting over a broad range of 
frequencies, transferring energy from the excited mode 
luq to other modes. The contribution to the value of 
the inverse of the quality factor, Q^^{ujo), of all these 
non-resonant TLSs will be proportional to A^^g_j.^^{u!o) = 

ElrJJ""' A(AS, Ag,^o) + E^:;r(.„)^(Ag,Ag,c.o) « 

2Pab\/i-Oco/i-^o : ^ quantity measuring the density of 
states which can be excited through Hi^t at frequency 
Wo, see Appendix IB] for details. 

For an excited mode populated with n^^Q phonons, 
Q~^{uJo) is given by Q~^(u;o) = /S.E/2ttEq, where Eq 
is the energy stored in the mode per unit volume, Eq ~ 
ric^ghuJo/twL, and AE is the energy fluctuations per cycle 
and unit volume. AE can be obtained from Fermi Golden 
Rule: 



27r ^ 27r 
— X nujo X —n^ 



fr2 

A- ° 



(5) 

and the inverse quality factor of the vibration follows. 
For finite temperatures the calculation of A^^_^^^{ujo,T) 
is done in Appendix |B] The result, valid for tempera- 
tures below the breakdown of the TLS approximate de- 
scription of the two- well potential (T ~ 5K— ), is dif- 
ferent from the one given in -i, because there the value 
of Q^^{iLjo,T) was interpreted as corresponding to the 
net energy loss of the studied mode (subtracting emission 
processes from absorption ones) , while in experiments the 
observed linewidth is due to the total amount of fluctua- 
tions, the addition of emission and absorption processes. 
So in this context dissipation means fluctuations, and 
not net loss of energy. The contribution of these kind of 
processes is thus 



[Q 



lOPot^/^w i — 

p 



1/4 2 



cotanh 
(6) 



T 



C. Contribution of biased TLSs to the linewidth: 
relaxation absorption 

This very general friction mechanism arises due to 
the phase delay between stress and imposed strain rate. 
In our context, for a given TLS the populations of its 
levels take a finite time to readjust when a perturba- 
tion changes the energy difference between its eigenstates 



"^^1^'^ . This time r corresponds to the inverse linewidth 
and is given, for not too strong perturbations, by the 
Fermi Golden Rule result 



r(£,T) 



/ecoth[e/2T] 



(7) 



where e = \J (Ag)^ -I- (A^)^ is the energy difference be- 
tween the levels. Notice that is not just the bare 
Aq appearing in the Hamiltonian of the system, but 
the net bias including the modification due to the cou- 
pling to the vibrational modes, A^ = Aq + £,kduk (Cfc 
is the corresponding coupling constant with the proper 
dimensions and duk a component k of the deforma- 
tion gradient matrix, defined recalling eq.Q as duk ~ 
{k,nk\{k/y/ujk){al. + ak)\k,nk), associated to a vibra- 
tional mode rife)). Eq.([7|) is valid in the range of appli- 
cability of the TLS description of the two well potentials 
(T < 5K), and for values of at such that r(e,r) < e. 
Therefore the energy levels £1,2 = T^-\/(A§)^ + (A^)^ ~ 
=Fi-\/(Ag)2 + (Ag + S^kdukY depend on duk, and to first 
order the sensitivity of these energies to an applied strain 
is proportional to the bias 



(A§ =F ikduk) 



d{duk) 



Ai 



(8) 



with a response of the TLS (x (A^/e)^. In — a de- 
tailed derivation is given, and the imaginary part of the 
response, corresponding to is also oc t/(1 -f- cj^t^), 
which in terms of A(uj) is the lorentzian peak at a> 
of fig(I33). 

The mechanism is most effective when r ~ l/i^, where 
uj is the frequency of the vibrational mode; then, along 
a cycle of vibration, the following happens (see fig. ([3])): 
When the TLS is under no stress the populations, due to 
the delay r in their response, are being still adjusted as 
if the levels corresponded to a situation with maximum 
strain (and therefore of maximum energy separation be- 
tween them, cf. £1,2(1*^)), so that the lower level becomes 
overpopulated. As the strain is increased to its maximum 
value the populations are still adjusting as if the levels 
corresponded to a situation with minimum strain, thus 
overpopulating the upper energy level. Therefore in each 
cycle there is a net absorption of energy from the me- 
chanical energy pumped into the vibrational mode. 

If on the other hand t ^ l/ui the TLSs levels' popu- 
lations are frozen with respect to that fast perturbation, 
while in the opposite limit r <C 1/w the levels' popula- 
tions are always in instantaneous equilibrium with the 
variations of £1,2 and there is neither a net absorption of 
energy. 

For an ensemble of TLSs, the contribution to Q^^ is 

46 



Q7j(c^,T) = PolViET) J de / duy/l - u^u x 



u„ 
LOT 



cosh^ (£/2T) 1 + {lot)^ 



(9) 
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FIG. 4: (Color online) Schematic representation of the levels' 
population evolution of an ensemble of 5 identical TLSs with 
a response time r ~ where to is the frequency of the 

bending mode excited, whose evolution is also depicted in the 
lower part of the figure. The delay r plus the bias Aq give 
rise to the relaxational energy loss mechanism of the mode, 
see text. 



The derivation of eq. © only relies in the assumptions 
of an existence of well defined levels who need a finite 
time r to reach thermal equilibrium when a perturba- 
tion is applied, and the existence of bias |Ag| > 0. This 
implies that such a scheme is applicable also to our ID 
vibrations, but is valid only if the perturbation induced 
by the bath on the TLSs is weak, so that the energy 
levels are still well defined. Therefore we will limit the 
ensemble to underdamped TLSs for whom r(e,T) < e. 
In eq.® the factor cosh~^ {e/2T) imposes an effective 
cutoff e < T, so that in eq.© one can approximate 
coth[e/2r] 2T/e, resulting in r(£) l/\/e, and thus 
the underdamped TLSs will satisfy e > [SOaby^uJ^T]'^^^ . 
For T 3> [32aby/ujco\'^ , which is fulfilled for typical sizes 
and temperatures (see Appendix [Cl for details) 



20Po7Ml + ^)(l 



Here Vamorph ~ Vtot/W was assumed 



') \EJ 



/4VT 



(10) 



IV. COMPARISON BETWEEN 
CONTRIBUTIONS TO Q-\ RELAXATION 
PREVALENCE. 

It is useful to estimate the relative importance of the 
contributions to Q^^ coming from the last two mecha- 
nisms. For that sake, we particularize the comparison 
to the case of the fundamental flexural mode, which is 
the one usually excited and studied, of a doubly clamped 
beam, with frequency ujq « 6.5{E/ py^^t/L"^ (for a can- 
tilever these considerations hold, with only a slight mod- 
ification of the numerical prefactors; the conclusions are 
the same). The result is: 



'Off' 



300ti/2 (3 - 



E 1 



fund 



L 



{l + u){l~2v) p Ti/2 

(11) 

For a temperature T = IK, the result is as big as 10^ 
even for a favorable case t = Inm, L = 1/xm, E = 50GPa, 
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FIG. 5: (Color online) Example of fit of Q-^^{ujq,T), eq.JTSJ, 
to experimental data. The data correspond to a Si resonator 
vibrating at ujq = 12.028 MHz, see ref.— . Although the pre- 
dicted order of magnitude for the dissipation and the sublinear 
temperature dependence is observed, the experimental trend 
follows a weaker temperature dependence than the predic- 
tion r^/^. For comparison a linear temperature dependence 
is shown (dotted line). See the text for more details. 



V = 0.2, p = 3g/cm"^. So for any reasonable temperature 
and dimensions the dissipation is dominated by the re- 
laxation mechanism, so that the prediction of the limit 
that surfaces set on the quality factor of nanoresonators 
is, within this model for the surface, eq. (fTO|) : 



surface 



(c^o,T)«Q7^[(c^o,T) 



2.1/2 



UJo 



(12) 



For typical values L ^ 1/im, t,w ^ 0.1/im, 7 ~ 5 eV, 
PaVamorph/Vtot ^ lO^^'J/m^, and T in the range ImK- 
0.5K the estimate for Q^^face ~ 1'^"'' gives the observed 
order of magnitude in experiments like — , and also pre- 
dicts correctly a sublinear dependence, but with a higher 
exponent, 1/2 versus the experimental fit 0.36 in or 
0.32 in — , see fig. ([5]) for an example. 



V. EXTENSIONS TO OTHER DEVICES 
A. Cantilevers, nanopillars and torsional oscillators 

The extrapolation from doubly clamped beams to can- 
tilevers ~ and nanopillars — is immediate, the only 
difference between them being the allowed (k,ix'(k)) val- 
ues due to the different boundary conditions at the free 
end (and even this difference disappears as one consid- 
ers high frequency modes, where in both cases one has 
kn ~ (2n -I- l)7r/2L). All previous results apply, and one 
has just to take care in the expressions corresponding to 
the of the fundamental mode, where there is more 
difference between the frequencies of both cases, the can- 
tilever one being wg""* « {E / pY^'^t/ L"^ as compared to 



the doubly clamped case. 



clamped 
^0 



&.b{E/pfl^t/L^ 
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Effect of the flexural modes on the dissipation 
of torsional modes 



The contribution from the TLSs + subohmic bend- 
ing mode environment to the dissipation of a torsional 
mode of a given oscillator can be also estimated. We will 
study the easiest (and experimentally relevant case 
of a cantilever. For paddle and double paddle oscillators 
the geometry is more involved, modifying the moment of 
inertia and other quantities. When these changes are in- 
cluded, the analysis follows the same steps we will show. 

Relaxation absorption. We assume, based on the 
previous considerations on the predominant influence of 
the flexural modes on the TLSs dynamics, as compared 
with the influence of the other modes, that the life- 
time T = of the TLSs is given by eq.®. The 
change in the derivation of the expression for Q^^ comes 
in eq.®, where the couphng constant £^1°'^^ is differ- 
ent, which translates simply, in eq.([9]), into substituting 
7^ ^ 7torsj s-iid the corresponding prediction for Q~J^ 



and the inverse quality factor 



Q 



-V rr., 20Po7t^„,,7^ (1 + ^)(1 - 2^) ^ p ^ 1/4 



T 

(13) 

where now tuo is the frequency of the corresponding flexu- 
ral mode. The range of temperatures and sizes for which 
this result applies is the same as in the case of an excited 
bending mode. 

Dissipation of symmetric non-resonant TLSs 
The modified excitation spectrum of the TLS's ensem- 
ble, ^off_j.es('^0 7 r), remains the same, and the change 
happens in the matrix element of the transition proba- 
bility of a mode |fco,«o) appearing in eq.®, {XkQ / ^/uJq)'^ . 
The operator yielding the coupling of the bath to the tor- 
sional mode which causes its attenuation is the interac- 
tion term of the Hamiltonian, which for twisting modes 
is (see Appendix [Dl for the derivation of eas. lfT4)) - (fT5| ) 



Triors 




(4 + a,.), (14) 



where /i — E/[2{\ + v)] is a Lande coefficient, and 
C = is the torsional rigidity. Again, Q^^{ujj) — 

AE /2TrEo, where the energy Eq stored in a torsional 
mode (t)j{z,t) = Asin[(2j — \)-k z / {2L)\ sm{ix} jt) per unit 
volume is Eq — A^ui'jpit'^ + i(;^)/48 {z is the coordinate 
along the main axis of the rod). Expressing the ampli- 
tude A in terms of phonon number Uj , the energy stored 
in mode | kj , rij ) is 

Eo{k„n) ^ 1 3 it' + w'){2n, + 1) , (15) 

the energy fluctuations in a cycle of such a mode is 

'^■-n,^^°ff*__(^„T), (16) 



27r ^ 

AE — — X huji X 



27r j'^htOj 



(Q"')off-rcs(^0,T) 



icotanh 



Y 



(17) 



For sizes and temperatures as the ones used for previous 
estimates the relaxation contribution dominates dissipa- 
tion. 



VI. FREQUENCY SHIFT 

Once the quality factor is known the relative frequency 
shift can be obtained via a Kramers-Kronig relation 
(valid in the linear regime), because both are related to 
the imaginary and real part, respectively, of the acoustic 
susceptibility. First we will demonstrate this, and after- 
wards expressions for beam and cantilever will be derived 
and compared to experiments. 



A. Relation to the acoustic susceptibility 

In absence of sources of dissipation, the equation for 
the bending modes is given by ~{12p/t')d'X/dP = 
Ed'^X/dz'^. The generalization in presence of friction 
is 



12p d'X 



= {E + X) 



d^X 



(18) 



Where x is a complex-valued susceptibility. Inserting a 
solution of the form X[z^t) = A^i-i^x-uji) ^ where k is 
now a complex number, one gets the dispersion relation 
bj = y/t'^{E + x)/(12p)fc^. Now, assuming that the rela- 
tive shift and dissipation are small, implying Re(x) << E 
, Im(k) << Re(k), the following expressions for the fre- 
quency shift and inverse quality factor are obtained in 
terms of x- 



Q-i = Aw/o) = -Im(x)/E 
5lj/uj =Re(x)/2E 



(19) 



Therefore a Kramers-Kronig relation for the susceptibil- 
ity can be used to obtain the relative frequency shift: 



Slo 



1 

2^' 



-P 



where P means here the principal value of the integral. 

B. Expressions for the frequency shift 

Relaxation processes of biased, underdamped TLSs 
dominate the perturbations of the ideal response of the 
resonator, as we have already shown for the inverse 
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quality factor. For most of the frequency range, to > 
[SOaby/U^T]^/^, Q-^{u},T) « Ay/T/uj , with A defined 
by eg. (fTO|) . The associated predicted contribution to the 
frequency shift, using eq. ipO)) . is 



(21) 



— (w,r) « log 



For low temperatures, uj '3> [SOaby/uJ^T]'^/'^ , the 
negative shift grows towards zero as Suj/lij{uj,T) 
7riog[T2/3/w], reaching at some point a maximum 
value, and decreasing for high temperatures, uj < 
[30af,7^:^r]2/3, as 6uj/uj{uj,T) - Even though 

the prediction of a peak in Slo /lo{T) qualitatively matches 
the few experimental results currently available — it 
does not fit with them quantitatively. 



VII. APPLICABILITY AND FURTHER 
EXTENSIONS OF THE MODEL. DISCUSSION 

As mentioned, the predictions obtained within this 
theoretical framework do match qualitatively experimen- 
tal results in terms of observed orders of magnitude for 
Q~^{T), weak sublinear temperature dependence, and 
presence of a peak in the frequency shift temperature de- 
pendence. But quantitative fitting is still to be reached, 
while on the experimental side more experiments need to 
be done at low temperatures to confirm the, until now, 
few results 

Applicability. The several simplifications involved 
in the model put certain constraints, some of which are 
susceptible of improvement. We enumerate them first 
and discuss some of them afterwards: i) The probability 
distribution P(Aq, Aq), borrowed from amorphous bulk 
systems, may be different for the case of the resonator's 
surface, ii) The assumption of non-interacting TLSs, only 
coupled among them in an indirect way through their 
coupling to the vibrations, breaks down at low enough 
temperatures, where also the discreteness of the vibra- 
tional spectrum affects our predictions iii) When temper- 
atures rise above a certain value, high energy phonons 
with 3D character dominate dissipation, the two-state 
description of the degrees of freedom coupled to the vi- 
brations is not a good approximation, and thermoelastic 
losses begin to play an important role, iv) For strong driv- 
ing, anharmonic coupling among modes has to be con- 
sidered, and some steps in the derivation of the different 
mechanisms, which assumed small perturbations, must 
be modified. This will be the case of resonators driven 
to the nonlinear regime, where bistability and other phe- 
nomena take place. 

The solution to issue i) is intimately related to a bet- 
ter knowledge of the surface and the different physical 
processes taking place there. Recent studies try to shed 
some light on this question—, and from their results 
a more realistic P(Ag,AQ) could be derived, which re- 
mains for future work. Before that point, it is easier to 



wonder about the consequences of a dominant kind of 
dissipative process which corresponded to a set of TLSs 
with a well defined value of Aq and a narrow distribution 
of Aq's of width Ai, as was suggested for single-crystal 
silicon Following a Q^^(T) ^ VT behavior is 
obtained for low temperatures T < Ai if r(Ag, T) < ui, 
and a Q~^T) - l/VT if r(Ag,T) > w, while at high 
temperatures T > Ai a constant Q~^{T) ^ Qq is pre- 
dicted for both cases. These predictions do not match 
better with experiments than the results obtained with 
P(Ag, Aq) ~ Pq/Aq, so issue i) remains open. From the 
point of view of those who attribute the origin of the low 
energy TLSs to the long range interaction between local- 
ized defects, a possible source of change for P(Ag, Aq) 
may be the decrease in size of the resonator below the 
correlation length of this interaction. 

We try now to have a first estimate of the temperature 
for which interactions between TLSs cannot be ignored. 
Following the ideas presented in we will estimate the 
temperature T* at which the dephasing time Tint due to 
interactions is equal to the lifetime t{T) = r^^(T) de- 
fined in eq.([7]), for the TLSs that contribute most to dis- 
sipation, which are those with e — (Aq) 
Ag - T. For them t-^{T) = r(T) w 40^6 
The interactions between the TLSs are dipolar, described 
by i/int = Z^ij t^i,20-!cr|, with C/i,2 = 6i2/rf2, 6i2 ver- 
ifying (&12) « 0, (I&12I) = C/o « I'^IE i&il. From 
the point of view of a given TLS the interaction af- 
fects its bias, 




(Ag)j Y.i^ii'^]^ causing fluc- 
tuations of its phase <5ej(t) (where t here means time 
and not thickness), which have an associated Tint de- 
fined by fcj (Tint)Tint ^ 1. Thcsc fluctuations are caused 
by those TLSs which, within the time Tint, have under- 
gone a transition between their two eigenstates, affecting 
through the interaction i?int the value of the bias of our 
TLS. At a temperature T, the most fluctuating TLSs 
are those such that e = ^(Ag)2 -f (Ag)2 ~ Ag - T, 
and their density can be estimated, using P(Ag,Ag), 
as riT ~ PokT they will fluctuate with a characteristic 
time t(T) « [iOab^/uJ^Vrj^'^ , so for a time t < t{T) 
the amount of these TLSs that have made a transition 
is roughly n{t) « PokTt/T{T). For a dipolar interac- 
tion like the one described above, the average energy 
shift is related to n{t) by — 6e{t) « Uon{t). Sub- 
stituting it in the equation defining Tint, and impos- 
ing Tint(r*) = T(r*) gives the transition temperature 
T* fa [6aby/uJco/ {UqPo)]'^- For example, for a resonator 
like the silicon ones studied in L = 6/im, t = 0.2/im, 
w = 0.3/xm, the estimated onset of interactions is at 
T* « lOmK. 

An upper limit T\^igh of applicability of the model due 
to high energy 3D vibrational modes playing a significant 
role can be easily derived by imposing Tiiigh — (^'^n- The 
frequency oj'^^ corresponds to phonons with wavelength 
comparable to the thickness t of the sample (do not con- 
fuse with time), uj^^ = pjt. The condition is 
very weak, as the value for example for silicon resonators 
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Charges in the oscillating part of the device. 

We write the relative distance as R.(t) — f (i) = R-o — f o + 
(5R(t) — Sr{t) and expand V{r,t), whose time-dependent 
part is approximately 



V{r,t) 



(5R(t) - Sr(t) 



eo|Ro - rol 



(24) 



FIG. 6: (Color onIine)Sketch of the distribution of charges in 
the device. When the system oscillates, these charges induce 
time dependent potentials which create electron-hole pairs in 
the metallic layer deposited on top of the beam, absorbing 
part of the mechanical energy of the flexural mode. See text 
for details. 



reads Thigh ~ 400/t, with t given in nm and Thigh in K- 
At much lower temperatures the two-state description of 
the degrees of freedom coupled to the vibrations ceases to 
be realistic, with a high temperature cutoff in the case of 
the model applied to amorphous bulk systems of T ~ 5K. 



VIII. DISSIPATION IN A METALLIC 
CONDUCTOR 

Many of the current realizations of nanomechanical de- 
vices monitor the system by means of currents applied 
through metallic conductors attached to the oscillators. 
The vibrations of the device couple to the electrons in the 
metallic part. This coupling is useful in order to drive and 
measure the oscillations, but it can also be a source of 
dissipation. We will apply here the techniques described 
in (see also — ) in order to analyze the energy loss 

processes due to the excitations in the conductor. 

We assume that the leading perturbation acting on 
the electrons in the metal are offset charges randomly 
distributed throughout the device. A charge q at position 
R interacts with an energy 



Vir,t) 



eomt)-m\ 



(22) 



with another charge g at a position r inside the metal, 
see fig. [6]. As the bulk of the device is an insulator, this 
potential is only screened by a finite dielectric constant, 
eg. The oscillations of the system at frequency luq modu- 
late the relative distance |R(t) — r(t)|, leading to a time 
dependent potential acting on the electrons of the metal. 

The probability per unit time of absorbing a quantum 
of energy ujq can be written, using second order pertur- 
bation theory, as 



58.59. 



T = J drdr'dtdt'V{r,t)V{r\t')lmx[r-r' ,t-t']e''^'>^*-'''^ 

(23) 

where Imx[r — r', t — t'] is the imaginary part of the re- 
sponse function of the metal. 



For a flexural mode, we have that Rq — fo • [(5R(t) — 

Sr{t)] has turned into [Xq — xo] ■ [5X{t) — 6x{t)] ~ t ■ 
A ■ sm(ujt), where t = t\ns + ^mctai is the thickness of the 
beam and A is the amplitude of vibration of the mode. 
Thus the average estimate for this case for the correction 
ofV{r,t) is 



SVir,t) 



t ■ A - s\n{ujt) 



L3 



(25) 



L being the resonator's length. 

The integral over the region occupied by the metal in 
eg. ([23)1 can be written as an integral over Tq and fp. The 
dielectric constant of a dirty metal in the Random Phase 
Approximation is — : 



Imx(r-r',<-i') = 



dudq 
(2^ 



giq(f-r')gii 



(26) 

where e is the electronic charge, D = hvpl is the diffusion 
constant, v-p is the Fermi velocity, I is the mean free path, 
and v w (k-ptnictid)'^ / (Tt-vf) is the one dimensional density 
of states. 

Combining ea. (P5)) and eq. pS)) . and assuming that 
the position of the charge, Rq is in a generic point in- 
side the beam, and that the length scales are such that 

^' ^motai, t <^ L, we can obtain the leading dependence 
of r in eq.(I2Sl) on L: 



DvL 



t 



metal 



(27) 



where we also assume that \q\ — e. The energy ab- 
sorbed per cycle of oscillation and unit volume will be 
AE = {2n/u!o)hujoTph/t''L = 2nnTpi,/t''L, and the in- 
verse quality factor Qph (t^o) will correspond to 



n-if ^ 1 ^rph 



1 



twL \pujlA^ 



(28) 



where Eq is the elastic energy stored in the vibration and 
A the amplitude of vibration. Substituting the result for 
r one obtains 



2n 



t 



(29) 



In a narrow metallic wire of width ^mctah we expect that 
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Typical values for the parameters in eq. ((29)) are fcp ^ « 
lA, A w lA, Z - ^mctai « lOnm ~ 10^ A, t « lOOnm 
w lO^Aand L w Ifim w lO'^A. Hence, each charge in 
the device gives a contribution to of order 10~^". 
The effect of all charges is obtained by summing over all 
charges in the beam. If their density is Uq, we obtain: 



2hn„ 



(30) 



For reasonable values of the density of charges, Ug = 
Iq^Jq ^ lOnm, this contribution is negligible, Q^^ < 
10-16. 

Charges in the substrate surrounding the de- 
vice. Many resonators, however, are suspended, at dis- 
tances much smaller than L, over an insulating sub- 
strate, which can also contain unscreened charges. As 
the Coulomb potential induced by these charges is long 
range, the analysis described above can be applied to 
all charges within a distance of order L from the beam. 
Moreover, the motion of these charges is not correlated 
with the vibrations of the beam, so that now the value 



of (5R(<) — Sr{t) has to be replaced by: 



(5R(t) - Sr{t) 



Ae 



(31) 



and the value of |R,o — fo| ~ L. Assuming, as before, a 
density of charges = the effect of all charges in the 
substrate leads to: 



2hL 



l{k 



F^mctal/ •^q 



0.3 



/(fcpi 



metal } 



where the second result corresponds to the fundamen- 
tal mode, ujo « 6.5(i/L^)-\/i?/p. For values L sa 1/xm, 
A « lA, fcpi « lA, I ^ imctai ~ lOnm and Iq ^ lOnm, 
we obtain Q^^ ~ 10^^. Thus, given the values of 
reached experimentally until now this mechanism can be 
disregarded, although it sets a limit to Q^^ at the lowest 
temperatures. It also has to be noted that this estimate 
neglects cancelation effects between charges of opposite 
signs. 



IX. CONCLUSIONS 

Disorder and configurational rearrangements of atoms 
and adsorbed impurities at surfaces of nanoresonators 
dominate dissipation of their vibrational eigenmodes at 
low temperatures. We have given a theoretical frame- 
work to describe in a unified way these processes, im- 
proving and extending previous ideas J*. Based on the 
good description of low temperature properties of disor- 
dered bulk insulators provided by the Standard Tunnel- 
ing Model 43,44,46 ^ g^^^ particular of acoustic phonon 
attenuation in such systems, we adapt it to describe the 
damping of ID flexural and torsional modes of NEMS as- 
sociated to the amorphous-like nature of their surfaces. 



Correcting some aspects of =, we have calculated the 
damping of the modes by the presence of an ensemble of 
independent Two-Level Systems (TLSs) coupled to the 
local deformation gradient field diUj created by vibra- 
tions. The different dissipation channels to which this 
ensemble gives rise have been described, focussing the 
attention on the two most important: relaxation dy- 
namics of biased TLSs and dissipation due to symmet- 
ric non-resonant TLSs. The first one is caused by the 
finite time it takes for the TLSs to readjust their equi- 
librium populations when their bias Ag is modified by 
local strains, with biased TLSs playing the main role, 
as this effect is oc [Ag/V(Ag)2 + (Ag)2]2. in terms of 
the excitation spectrum of the TLSs, it corresponds to a 
lorentzian peak around lu = 0. The second effect is due 
to the modified absorption spectrum of the TLSs caused 
by their coupling to all the vibrations, specially the flex- 
ural modes, whose high density of states at low energies 
leads to subohmic damping —i^i^. A broad incoherent 
spectral strength is generated, enabling the " dressed- by- 
the-modes" TLSs to absorb energy of an excited mode 
and deliver it to the rest of the modes when they decay. 

We have given analytical expressions for the contri- 
butions of these mechanisms to the linewidth of fiexu- 
ral (eqs.dni), pI7|) ) or torsional modes fegs. fT^ . (fT7|) ) in 
terms of the inverse quality factor Q~^{ujo) = Awq/wo, 
showing the dependencies on the dimensions, tempera- 
ture and other relevant parameters characterizing the de- 
vice. We have compared the two mechanisms, concluding 
that relaxation dominates dissipation, with a predicted 
Q^^{ujo, T) ^ T^l"^ jujQ . Expressions have been provided 
for damping of fiexural modes in cantilevers and doubly- 
clamped beams, as well as for damping of their torsional 
modes. 

Analytical predictions for associated frequency shifts 
have been also calculated fea. (|2T|) ). Some important suc- 
cesses have been achieved, like the qualitative agreement 
with a sublinear temperature dependence of Q~^{T\ the 
presence of a peak in the frequency shift temperature de- 
pendence 5ui/u:{T), or the observed order of magnitude 
of Q~^{T) in the existing experiments studying flexural 
phonon attenuation at low temperatures ^iS^. Neverthe- 
less, the lack of full quantitative agreement has led to a 
discussion on the assumptions of the model, its links with 
the physical processes occurring at the surfaces of NEMS, 
its range of applicability and improvements to reach the 
desired quantitative fit. 

Finally, we have also considered the contributions to 
the dissipation due to the presence of metallic electrodes 
deposited on top of the resonators, which can couple to 
the electrostatic potential induced by random charges. 
We have shown that the coupling to charges within the 
vibrating parts does not contribute appreciably to the 
dissipation. Coupling to charges in the substrate, al- 
though more significant, still leads to small dissipation 
effects fea. (|5^ ). imposing a limit at low temperatures 
~ 10"^, very small compared to the values reached 
in current experiments ^i^^. 
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ments of the rod as a function of the position along its 
length, z) obeying — 



EL 



EI. 



—ptw 
—ptw 



(Al) 



APPENDIX A: CALCULATION OF THE 
SPECTRAL FUNCTION FOR THE BENDING 
MODES 

The starting point is the Hamiltonian of ref. ^, H = 
^oCk +l<^zdiUj, where diUj is a component of the defor- 
mation gradient matrix. In the case of the bending modes 
of a rod of dimensions L, t and w, and mass density p, 
there are two variables X{z),Y{z) (transversal displace- 



J 



(where ly = t^w/12, and I. ~ w^i/12), so that there 
are plane waves X{z,t),Y{z,t) ^ gi(kz-ujt) ^ with a 



quadratic dispersion relation, 0Jj{k) 



X k . One 



can thus express X{z)^Y{z) in terms of bosonic opera- 
tors, for example X{Q) = Y.k \/ 2ptwLu, (4 + We 
can relate this variables to the strain field diUj through 
the free energy F: 



o2 x'' 2 /^^^-^ 2 1/" /"I 

F.., = ^ / dzEIy[^) ^^J dzj dS-XY^ul + t^Y.^^'k dz I dSi^X + 9^)ul (A2) 

i i,k 



extracting an average equivalence for one component , term in the Hamiltonian is then 
-Uy K. 2^jEIy/{i\ + 18n)twd'^X/dz'^. The interaction 



k V^'= 



ha. 



ij k 



Eh 



(3A + 18^)tw V 2ptwhL 



(A3) 



So we have approximately 9 times the same Hamiltonian, 
once for each Uij , and the corresponding spectral function 

I 



J{uj) will be nine times the one calculated for 



Hint = TicTz ^ \—={a\ + ak) ~ Tia, 




EL, 



(3A + 18^i)tw V 2ptwhL 



/UJk 



(al+ak) 



(A4) 



For a Hamiltonian of the class H = Aocr^: -I- by J(w) = '^J^k^t^i'^ ~ ^k), so that in our case the 
huz X]fc ^kicil. + o-k) the spectral function J{uj) is given expression for it is 



EL, 



1 



(3A + 18ii)tw V 2ptwhL y/Uk 



S{uj - UJk) 



(A5) 
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Taking the continuum limit (-^ J^k ^ ^ / '^^)' 



J{lo) 



2L 
(2^ 



dk 



27 



1 



Ag / 

Ao y (3A + 18fi)tw V 2ptwnL 

r 



■(5(w - ujk) 



(A6) 



Using the dispersion relation ojj (k) = y xk^ — cxk^ 
we express the integral in terms of the frequency: 



Jico) 




AX 



Eh 



Ao Y (3A + I8n)tw V 2ptwhL 



— 6[uJ-UJk) = jTT-yl 




Eh 



r 



(3A + I8fi)tw V 2ptwhL 
m 



'/flcx(^) is just 9 times this, eq.([2]). 



APPENDIX B: DISSIPATION FROM 
OFF-RESONANCE DRESSED TLSs 

We follow the method of ref. The form of A{lo)^ 
the spectral function of a single TLS, for frequencies 
w A and cj ^ A can be estimated using perturba- 
tion theory. Without the interaction, the ground state 
\s) of the TLS is the symmetric combination of the 
ground states of the two wells, and the excited state is 
the antisymmetric one, \a). We will use Fermi's Golden 
Rule applied to the subohmic spin-boson Hamiltonian, 



(4 



ak) 



'£k^^ik)aiak, 



where a^. is the annihilation operator of a bending mode 
k. Considering only the low energy modes uj^k) <C A, to 
first order the ground state and a state with energy uj{k) 
are given by 



Iff) 
\k) 



Xk-' 



2A 



4 1 a) 



4k> 



\k^ 
'2K 



(Bl) 



We estimate the behavior of A(io) by taking the matrix 
element of CTz between these two states, obtaining (re- 
member that u;(fc) oc fc^) 



the ground state \g) and an excited state \k) can be writ- 
ten as 

\k^/^ t 



(B3) 



The matrix element (OIctzI/c) is ~ A^-p-^A_^ leading to 

. , cj(fc)>A. (B4) 



1. Value of ^*?f^ 



Now we will add the contributions of all the 
non-resonant TLSs using the probability distribution 
P(A5, Ag) = Po/Ag For the case of weak cou- 

pling, CKfc < 1/2 and wq > {2ai,Yu!co, which is the one 
found in experiments, one has 



^ofF*(^0) 



S Ag (Ag)2 



dAg 



/\x 3 7/2 



obtaining the result A)^q_^^^{ijJo) « 2Pah\J '^col ^a- 



){k) < A. (B2) 



2. The off-resonance contribution for T > 



The expression in the numerator is proportional to the 
spectral function of the coupling, J(iS) = a\,^Joj^\/u) ■ 
Now we turn our attention to the case a;(fc) 3> A, where 



Using the same scheme, the modifications due to the 
temperature will appear in the density of states of absorp- 
tion and emission of energy corresponding to a " dressed" 
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TLS, A{A, ui, T). Now there will be a probability for the 
TLS to be initially in the excited antisymmetric state, 
|a), proportional to exp[— A/fcT], and to emit energy Ttw 
giving it to our externally excited mode, \k,n), thus com- 
pensating the absorption of energy corresponding to the 
opposite case (transition from |s)|fc,n) to |a)|fc,n— 1)), 
but contributing in an additive manner to the total 
amount of fluctuations, which are the ones defining the 
linewidth of the vibrational mode observed in experi- 
ments, fixing the value of (5~^(a;,T). The expression 
for j4(A, w, T) is given by 

A{^,T) = |(i|a,|/)|2e-|^^[fta; - {E^ - E^)] 

i f 

(B5) 



We consider a generic state \ia) = \a)\kini, ...,kjnj, ...) 
or \ig) = \s)\kini, ...jkjUj, ...), and states that differ 
from it in hojj, \fa±) = \a)\kini, ....kjUj ± 1,...) and 
\fs±)\s)\kini, ...,kjnj ± 1, ...). As for T = 0, we will cor- 
rect them to first order in the interaction Hamiltonian 
Hint = ^^cTz^k V'^i'A: + '^k), and then calculate the 
square of the matrix element of cr^, |(i|fT2|/)p. 



Elements \{ia,s\o'z\fa+,s+)\'^ correspond to absorption 
by the "dressed" TLS of an energy hu)j from the mode 
kj, while elements \{ia,s\(^z\fa-,s-)\'^ correspond to emis- 
sion and "feeding" of the mode with a phonon JkVj. The 
expressions for the initial states are 



\is) = \s)\kini,...,kjnj,...) \i 



\kini,...,kjnj,...) + ^ t;, . , oa \a)\-hni - I...) 



ki3ni>0 



TlUJi + 2A ' 



\ia) = \a)\kini,...,kjnj,...) \a)\kini, ...,kjnj, ...) + ^ 



ki3ni>0 



hwi - 2A 



\s}\...kini - 1...) 



Vfei 



fkJi + 2A 



(B6) 



and, for example, the state |/a+) is given by 



\fa+} = \a}\kini,...,kjnj + 1,...) \a)\kini, ...,kjnj + 1, ...) + ^ 



hwi - 2A 



\s)\...kini - L.-kjUj + 1...) 



huj - 2A 



hX^/ {rij + 2)ujj 
huj-i + 2A 



|s)|...fc,-n,-+2...) 

I 



hui + 2A 



(B7) 



with similar expression for the rest of states. The value 

of \{ia,s\(Jz\fa+,8+)\ (absorption) is \{ia.sWz\fa+,s+)\ = 



A-^(nj-|-l)wj4A 
(?i(^j)^-4A^ 



with : 



0, 1, while for emission the 



with rij = 1, .... 



result is \ {ia.s\<^z\.fa-,s-)\ 

Taking the limits we considered at T = {hiVj <C A and 
Uwj A) the results are the same as for T = for ab- 
sorption, except for a factor {nj + 1), and we also have 
now the possibility of emission, with the same matrix 



element but with the factor n,-: 



Absorption 
Emission 



(nj + l)ai,^t^eoA 



FZJTT^ , nuj > A 



~ « ^ (B8) 

Now we have to sum over all initial states, and noting 
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that the first order correction to the energy of any eigen- 
state is 0, the partition function, Z, is easy to calculate, 



J 



everything factorizes, and the result is, for example in 
the case huJi <^ A: 



OO OO t; / , 



A2 



{uj + l)e ST 



A2 1 



i^j^m—O rij — l 



A2 



A2 



1 — e fcT 



r 



(B9) 



In v4abs(A, Wj, r) we have added the contribu- 
tions from the matrix elements | (is |crz|/s+) pe'feT + 
\{iaWz\fa+)\^e-kT = |(is|crz|/s+)P[e^ +e^] . The sum 
of exponentials cancels with the partition function of 
the TLS (appearing as a factor in the total Z), leading 



J 



in this way to the expression above (the same applies 
for Aem(A, — Wj, r)). The total fluctuations will be 
proportional to their sum, which thus turns out to be at 
this level of approximation cx cotanh[ha;j/kT]: 



In fact, this result applies for any other type of modes, 
independently of its dispersion relation, provided the cou- 
pling Hamiltonian is linear in CTz and (alj, + a,k) and ev- 
erything is treated at this level of perturbation theory. 
Moreover, it can be proven that if one has an externally 
excited mode with an average population (rij) , and fluc- 
tuations around that value are thermal-like, with a prob- 



ability (X exp[ 



j)\/kT)], one recovers again the 



same temperature dependence, oc cotanh[hwj/kT]. 





O)cotanh 






TCA) 


CO 






1 


2 3 



(BIO) 



A/(D 



APPENDIX C: DERIVATION OF Q-[(LJo,r), 
EQ.dlOD 

As discussed after eq.®, we have to sum over under- 
damped TLSs, e > [SOaby'ujcoT]'^^'^ , using the approx- 
imation for r, r(e,T) ~ SOaby^LUcoT/ \/e. Moreover, if 
^0 > r(e = [30aby^T]'^/^,T) then in the whole in- 
tegration range luq » T{e,T) ^ ujQT{e,T) ^ 1 (see 
fig-Q, so that loqt/[1 + (wqt)^] ~ l/(wor), see fig.©. 
Q-f(wo,r) follows: 



ET 



de 



du X 



VI - r(e, T) 



Wo 



(CI) 



For temperatures T ^ [SOaby/uJ^T]'^^^ , which holds 
for reasonable T and sizes, the integral, which renders a 



FIG. 7: (Color online) Evolution with of the different 
quantities determining the approximations to be taken in the 
integrand ujot/[1 + (loot)'^] of eq.®. 



result of the kind Q-^1{luo,T) « {1/ujo){AVT - BT^^^), 
can be approximated by just the first term, obtaining 
eg. pU)) . In any case for completeness we give the expres- 
sion for B: 



B 



500Po7i^/3 (1 + - 2t/)4/3 / p 

£;ii/3(3_ 5i.) [e 



1/3 



(C2) 



Also for completeness we give the result for higher tem- 
peratures, although for current sizes the condition T(s — 
[SOab^yuJcoT]'^^^ , T) > implies values of T above the 
range of applicability of the Standard Tunneling Model. 
Now for some range of energies F(e,T) > wq, and the 
range of integration is divided into two regions, one where 
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luqt ^ 1 and one where the opposite holds: 



ET 



du 



[30abV^r]2/3 



dE uJoTie,T)+ [ 

J[16q. 



dE- 



The final result is Q~]{uja,T) w -~7Poj^ujo/T + 
AVT/ujo - CT/loI, with A defined by eq.HU]) and C by: 



1500Po7^ {l + vY[l-2yf 



1/2 



(C4) 



All the results for Q^.^{ujo, T) have to be multiplied by the 
fraction of volume of the resonator presenting amorphous 
features, VamorphI {twL). 



APPENDIX D: DERIVATION OF EQS. p4l) -(fT5 l) 



To derive the interaction Hamiltonian, eq. (jl4p . note 
that in terms of the deformation gradient matrix diUj 
it must be Hn^t — "fazdiUj, where the deformations are 
caused in this case by twisting of the resonator about 
its main axis. The twisting modes correspond — to the 
rotation angle around the longest main axis (j) obeying 



the wave equation 



(Dl) 



so the variable cj) can be expressed in terms of boson 
operators 



\z=0 



2ptwL 



+ (D2) 



/LUk 



To obtain in terms of the modes' operators Uk an approx- 
imate expression for diUj, we relate du to dcf) through the 
expressions for the free energy of the rod in terms of both 
variables, 



dz / dSAfi 



dz 



duy 
dz 



— I dz I dSSfi 



/ dux\' 



(D3) 



r 



and the approximate relation du |2=o= dux/dz = 
^ C I {8fitw)d(j) \z=o is found. This relation together with 



eq. (|D2p lead to the stated result, eq. lM]) . 

Calculation of Eq.To calculate classically the en- 
ergy stored of a torsional mode (j)j{z,t) = Asin[(2j — 



J 



l)TTz/{2L)]sin{LUjt) we just calculate the kinetic energy 
in a moment where the elastic energy is zero, for exam- 
ple at time t = 0. If an element of mass is originally 
at position (x,y,z) {x,y transversal coordinates), with a 
torsion 4>(z,t) it moves to 



+ ■ cos 



y^ ■ sm 



yx^^+y^ 



, z 



(D4) 



The kinetic energy at time t = is 
Eo = dz dx I dy ■ -p 

Jo J-t/2 J-w/2 ^ 



dt 



(D5) 



Substituting the expression for r{t) in the integrand, one 
arrives at 
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-En 



dz 



t/2 



dx 



t/2 



w/2 ^ 

dy ■ -p^^wj sin^ 

-w/2 ^ 



(2j - l)7r 



2L 



(D6) 



In terms of the creation and annihilation operators its amphtude is = ^2/2 = ?i(2n + l)/[2Lp/a;j]. Sub- 
(j)j{z,t) = 2Lpii^ (a]+a.,-)e''^^J^^"*\ so the mean square of stituting this in eq. (jD6|) the eq.lUS]) for Eq is obtained. 
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